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ABSTRACT 


The report investigates the viscous flow between two parallel disks rotating in 
the same direction with the same velocity. The fluid enters the space between the 
two disks at a certain radius in the radial direction, Secause of the shear forces, 
it assumes a rotating motion with about the velocity of the disks, The centrifugal 
forces then bu‘ld up a pressure increase in the radial direction, The arrangement 
corresponds to a centrifugal fluid pump, which may be advantageous if cavitation 
is a problem, 

The general equations of viscous flow are simplified by the assumption that 
the pressure difference normal to the disks is negligible (boundary iayer assump- 
tions), One obtains a system of parabolic partial differential equations. For large 
radii the deviation from rigid body rotation (with the angular velocity of the disks) 
is smail. The iin: ...7s i equations which then result are s. - analytically, The 
velocity profiles depend pon a parameter containing the kinematic viscosity, the 
angular velocity and the distance of the disks, but not the radius. 

The non-linearized parabolic differential equations are approximated by a 
difference acheme and solved numerically. The results are given in non-dimen- 
sional form with the entrance velocity and the distance of the disks as parameters, 


Furthermore, the efficiency of the pump is computed from the gain of the total 


pressure and the torque at the shaft of the rotating disks. 
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I. INTRODUCTION 


Investigated in this report is the viscous flow between two parallel disks, 
which rotate with the same constant angular velocity in the same direction. The 
fluid enters the space between the diaks in radial direction through a cylindrical 
surface (Figure 1), It is set into a spiral motion whose tangential and radial 
velocity components are due to the friction forces and to the centrifugal forces 
respectively. The pressure will increase with the radius because of the centrif- 
uga! forces, thus the arrangement can serve as a centrifugal pump. One might 
suspect that the absence of blades leads to low values of the efficiency, but if the 
blades are spaced very closely and the mass flow is low, the efficiency may be 
quite acceptable, On the other hand, such a pump would not encounter the problem 
of cavitation and also its characteristics over a wide operatin’ -ange might be 
more favorable tha: in a2 conventional turbo pump. These advantages were recog- 
nised by Mr. S.!. Hasinyer, Mr, L.G, Kehrt and Dr. J. P. von Ohain of the 
Thermomechanics Research Branch of this laboratory, Technical details and 
experimental results will be published by them in a future report. 

The present mathematical analysis will determine the velocity distribution 
between the disks, the pressures, the torque applied at the shaft and the resulting 
efficiency. 

There exists an extensive literature on the single disk rotating in an infinite 


medium*, Assuming that the disk extends from radius zero to infinity, it is 


* Advances in Applied Mathermatics, Academic Press, Inc.. New York, Vo! IV, 


1956, p. 166. F.K, Moore: Three-Dimensional Boundary Layer Theory. 


possible to reduce tne partial differential equations governing this problem to a 
sy3tem of ordine zy differential equations by means of a similarity hypothesis, But 
an approach of this kind is not feasible in the present case because of the presence 
of the seccid disk and the fact that the flow enters the space between the disks at a 
radius different from zero. However, another simplification can be made here. 
For any practical application the distance between the disks must be very small, 
for otherwise the friction would not be sufficiently effective to produce a tangential 
flaw, But then the approximations made in boundary layer theory apply reducing 
the original system of Navier-Stokes equations, vhich is elliptic, to a system of 
parabolic differential equations. They are derived in Section II. 

A first insight into the basic characteristics of the proble 1 is obtained if one 
assumes that the velocities of the fluid relative to the disks are small so that 
quadratic terms >: cogligible, This is a condition whic: ‘ell satisfied for 
large radii; tt 3 idea is carried out in section III, It gives valuable insight into 
the general character o/ the velocity profiles; these profiles are governed by one 
dimensionless pararaeter which contains the angular velocity, the vi=cosity and 
the distance, but not the mass flow and the radius, Naturally this parameter is 
important for the practical design. The solutions thus obtained can be considered 
as an asymptotic expression for the solution of the actual non-linear system. 

In the - icinity of the entrance radius the linearization is not justified, in par- 
ticular the tangential velocity relative to the disks is certainly not small. This 
region is most critical from the point of view of cavitation, and also the principal 
losses are encountered there, Therefore, the system of differential equations has 


been integrated numerically. The methods and the results are shown in section 


IV and V and in an appendix. 


Nl, STATEMENT OF THE PROBLEM 


We use a systcin of cylindrical coordinates r and z (Figure 1), where the 
z axis coincides with the axis of symmetry and r is the distance from it. We 


denote by 


u_- the component of the velocity in the radial direction 


v_ the component of the velocity in the tangential direction 


measured in a system of coordinates fixed in space 

w the component of the velocity in the axial direction 
2d_ the distance between the two disks 

ry the radius of the entrance hole 

p the density 

p the viscosity, v = the kinematic viscos?’ * 

p th. pre sure 

w she cons.ant angular velocity. 


An incompressible axi-symmetric steady viscous flow is described by the Navier- 


Stokes equations 


2 2 

8 ° 
um twee yey ee ey gee 
dr az 2 or r 2 
az or 


By applying the boundary layer approximations, this system of equations is reduced 


to 
8u Vv du lL ad 9¢ 
une = + = <P + gu 
or r oz p adr 822 
(2. 3) 
, 2 
u = ou tw oy z v oy 
dr r dz 
' 8z 
8u u bw 
ama Uy OX ra 
at 5, 70 (2, 4) 


It is convenient to introduce instead of v, which is the tangential velocity in a 


fixed system of coordinates, the tangential velocity V relative to the disks 


—— ie 


Vaverw (2, 5) 
One then obtains 
2 22 2 
pet, Mob evrete es py S85) ae yy ee (2. 6) 
or r z Pp dr 82° 
2 
ae eee See ve $y eee (2.7) 
or r az 922 


Along the disks, given by plane z = +d, one has the boundary ccaditions 


u (r, +d) =0 
(2, 8) 


V(r, edj)=0 


To determine the boundary value problem comrletely, one must also prescribe 


the velocity dist: ioution in the entrance cross section. One has 


u(r,, z)s Uy 


(2.9) 


V(r.) ge rw 


Ill, LINEARIZED TREATMENT AND ASYMPTOTIC SOLUTION 


Omitting quadratic terms in u, V and w one obtains from Eqs (2. 6) and (2. 7) 


9° i 4 
¥v OM g2yve —— SZ. rw (3. 1) 
W 92¢ pw dr 
2 
vy bv 
u+55 a2 (3. 2) 


dr 


assumption tha: the pre«sure does not depend upon z. The unknown w does not 


In these equations i oR does not depend upon z becauc che boundary layer 


appear, and so it is possible to solve these equations separately; w is determined 
from Eq (2.4) afterwards, It is scmarkable that :n these equations -uere appear 
no par.‘ai derivatives with respect to ». This means that, except for the expres- 


sion 


— SP. rw = F(r) (3. 3) 


which depends upon r and ia undetermined so far, the equations can te solved 
for each value of r independently. The expression (3, 3) is the only innomogeneous 


term that occurs in the equations and in the boundary conditions; therefore, it will 


enter the solution for u and V asa factor which depends only on r. This factor 


is determined by the condition that the mass flow between the disks is constant for 


everv cross section r=const. Thus u and V as wellas ; 4 can be found, 


After elimination of u one obtairs from Eqs (3.1) and (3. 2) 


a fh 
JY AV gave Flr) (4,4) 
29% dz 


Setting \= = u one has finally 


av 4 
q + 4e 


V = 2a? F(r) (3.5) 
dz 


As the most general solution of Eq (3.5) that satisfies the symmetry conditions at 


the plane 270, onc C12}. 
V=A, sinh (az) sin (az)+A, cosh(az)cos(az)+ +t) (3. 6) 


where 4, and A, are two constants which will be determined presently, The 


boundary conditions (2,8) together with Eq (3.2) give for 2 =2d 


(3. 7) 


Hence, by inserting Eq (3. 6) 


with 


_— 2 sir’. ,ad)- sin (ad) 


1 cosh(Za0) + cosf2ad) 


s 2 cosh(ad)+ cos (ad) 
2 cosh (2ad) +cos (2cd) 


Figure 2 shows k, and k, as a function of ad 


The shear force at the wall is proportional to dV/dz, One finds 


dV 
dz 


Let the flow of liquid (measured by volume) be given by Q. 


a surface r=const extending between the disks is given by 


+d td 2 
° Vv ° Vv 
Qz=2nr f udz= === | —_ dz = 25.2 
od @ od dz a 


Hence with Eq (3. 10) 


Q=-~- TS. F(r)eg (ad) 


—zsea Flr) : (k “k.) cosh(az) nin (az) +(k, +k) sinh (as) cos (aa) | 


(3. 8) 


(3.9) 


(3, 10) | 


The flow through 


+d 


ed 


(3.11) 


(3. 12} 


where: 


a) sinh (2ad) - sin (2d) 
glad) = —osh(zad) #cos (zed) (3. 13} 


Figure 3 shows g(ad) as a function of ad 


From Eq (3. 13) we find the function F(r) 


aQ 2q 
F(r)a¢ wereglad) © (3.14) 
where 
aQ 
4* 25. glad) het?) 
| 
One thus obtain: finally .or the velocity components, expressed in terms of the 
mass flow 
| Vs 4 [*, sinh(az) sin ( z)+k, cosh(az) cos (az) 1] (3. 16) 
| (3. 16) 
| us afr, cosh(az) cos (a:, = k, sinh (az) sin (<2)] 
Inserting these expressions into the continuity equation Eq (2.4) one finds w=0. 
The velocities in the middle plane between the two disks are found by setting z= 0. 


(3, 17) 


Neco [x -1] 


Figure 4 shows the radial velocity component for different values of ad as a func- 
tion of 2/d for q/r =1. For small values of ad, the profiles have a near para- 
bolic shape with maximum velocity in the middle (z=0). For ad = 1/2 the profile 
is flat in the riddle, for higher value of ad the flow is mor2 and more concentrated 
in the neighborhood of the disks. 

The tangential velocity V is represented in Figure 4 for q/r = 1. For small 
values of ad it also has a parabolic shape, but with increasing ad the slope in the 
vici::ity of the disks becomes steeper and the relative velocity in the middle plane 
ayaroaches -1. 

Actually the profiles for high values of cd must be viewed with considerable 
caution, One sees that in this case the velocity u may become n xgative in the 
middle; i.e. for such cases one has & radial outflow close to the disks with a 
radial inflow in the middle. Obviously it would take a rathe. tial arrangement 
to produce such: fiow pattern physically, These reaults ought to be disregarded 
for technical applications. Important in any casc is the role played by the dimen- 
sionless parameter ad which determines the character of the profiles. This 
character is neither influenced by the mass flow Q nor by the radius r. Fora 
given w anda given v, i.e., a given fluid, the only quantity that can be influenced 
in this parameter is the distance of the disks. If one wants to use high angular 
velocities, one must try to lower the value of d by skillful design. 

The pressure distribution is obtained from the function F(r) in Eq (3.5). 
Inserting Eq (3. 14) into (3. 3). 


one finds 


SP = pw [ro- “9] (3.18) 


and by integration 


2 
p= ra [ oa - 2qlnr + const ] (3. 19) 


The momentum required at the shaft of the rotor can be determined in two 
different ways, either from the moment of the shearing forces on the inner side 
of two disks, or by applying the law of conservation for the moment of momentum, 


One has in the first case 


r 
M =2-2n-y f 2” (SY ) a (3. 20) 
r dz 
i zed 
in the second 
2 td 
Mu2ar'p f uevdz (3.21) 
ed 


For the numerical evaluation the second one is preferable, 


IV. NUMERICAL SOLUTION OF THE SYSTEM 
OF EQUATIONS (THE INLET PROBLEM) 

The approximation shown in the previous section is unsatisfactory in the 
vicinity of the entrance, where the velocity components are certainly not small, 
Moreover, the linearization causes all r-derivatives to vanish, i.e. it introduces 
such *% strong change of the character of the flow equations, that it becomes 
impossible to satisfy the boundary conditions at the inlet. But the inlet is im- 


portant for it accounts for most of the losses. 


10 


As was mentioned above, the introduction of the boundary lz,er assumption 


gives us a parabolic system of differential equations where lines r = const are the 


characteristics, therefore, it is possible to proceed in the numericai solutions 


from one line r=const tu the next one without any need to go hack, 


Thus, we may 


adopt a numerical scheme for the integration of the parabolic differential equations. 


For our computations we assumed that the rotating disks lie in horizontal 


planes, and introduce a new system of coordinates x, y, whose origin lies in the 


entrance cross section at the lower disk. Thus the x-axis coincides with the 


lower disk, We have 


We introduce dimensionless variables in the system of “ferential equations 


(2.6), (2.7), and (2. 4) 


r. Yr. 

1 1 
Veve y= 2 qa- o 
- Ww 
WwW % 


11 


One thus obtains 


- oa - ¥ a 2. 
oe = ae 
“ag > Vrex $204 6 3 © Fea 2 


av Vv V z= 
[~+—pt2ajevn VV, Ov (4, 2) 


where for abbreviation 


F(r) a re —L— SP (4, 3) 
dr 


The functions G ind Y¥ and W are given at the entrance (x=0)}. The channel 
extends from y=0 to y¥= 2d. The boundary condition at y = 2d can be replaced 


by a symmetry coudition at y = d. Thus, one has as boundary conditions 


g(z,0)=0 alka) 9 
ay 

Viz,0-0 2¥bed L5 (4, 4) 
dy 


w (x, 0) =0 w(x,d) = 0 


The natural method of solving such a system numerically is the introduction 


of a difference scheme. Here caution is needed from the point of view of stability. 


12 


a 


The coefficients of the derivatives with respect to x, in Eq (4,2) may become very 
sma)l for large values of x ana it also is small in the vicinity of y= d, but then it 
is necessary to use a very fine mesh size in the x direction, if one applies a 
direct" differance scheme where the derivatives with respect to y are computed 
along the line x = const which is already known, Otherwise the procedure would 
become unstable, 

For this reason it was decided to use an "inverse'' difference procedure where 
the lerivatives are formed along the line x = const for which the state is to be 
computed, This brings about a rather severe complication, for the stat« along 
the new line x = const is determined by a system of simu!taneous equations, 

Since the original differential equations are non-linear, this system of equations 
is also non-linear, At first the non-linear system was solved directly by an itera- 
tion method, Since this method proved to be inconvenient, >2proach was modi- 
fied in the follow ng manner: In going from one line X = const to the next one we 
computed the changes of 41, V and W rather than the values themselves. If a 
difference procedure is admissible at all, these changes ought to be small and 
second order terms in these changes are negligible. In other words we use a 
linearization which considers as a basic approximation the values at the line 

x = const which has Leen computed previously. 

The state along the new line can then be obtained by a suitable super-position 
of particular solutions of the linear equatioa which arises in this manner, The 
linear system of equation which arises in this process has a rather simple struc- 
ture insofar as the matrix of systems has elements that are different from zero 


on only a few lines that are parallel to the main diagonal. This brings about 


certain simplifications in the inversion of the matrix, For further details of the 
13 


numerical analysis see the appendix. 
V. RESULTS AND CONCLUSIONS 


The numerical procedure described above gives a nearly complete description 
of the flow pattern. Assuming constant velocity Uy at the entrance, we give the 
results for a set of four initial radial velocities is 1,0, 0.5, 0.25and0.1. The 
initial tangential velocity Yi is always ve =-1 because of Eq (2.5). For the 
dimensionless distance between the disks we chouse d=0.5, 1.0, 1.5, and 2. 

Figures 5 thru 20 show the dimensionless velocity profiles versus z/d with 
r/t; as parameter, These pzofiles have a common characteristic: at the wall the 
inlet velocities Uy and Yo which are constant over the cross section, are immedi- 
ately reduced to zero because of the boundary conditions. The reduction of the 
radial velocities at the wall causes the velocities in the middl. :o increase. There- 
fore the radial velr ..ty '. the profiles close to the entrance c..us section over- 
shoots the entra: -e profi.© inthe middle. The effect decreases with increasing r 
because of the increase of the available cross section with r, finally the profile 
approaches the form given by the linearized theory, Figures 21 thru 24 show this 
behavior again; here the radial velocities ia in the median plane of the disks 
(z=0) are plotted versus r Ir, . The average radial velocities are obviously de- 
termined by the condition of continuous flow. 

The average tangential velocity is directly connected with the moment M of 
the shear forces and thus with the work that must be performed to drive the pump. 
Further.nore it gives the main contribution to the dynamic pressure of the fluid 


particles as they leave the pump, thus it influences strongly the efficiency. As 


mentioned before shear forces can be obtained from the slope of the velocity 


14 


profiles at the wall Eq (3.18); but it appears to be preferable (from the numerical 
point of view) to compute it from the radial and tangential velocities Eq (3.19). 

The pressure distribution is computed from the function F(z) which is found 
in the numerical procedure as a function of r/r, . Here an integration in the 
direction of r is required, In our figures the contribution to p due to the term r 
in Eq (4. 3) is shown separately, (as a straight line), The pressures are given by 
the difference of the curve r(#) from this straight line. 

Properly speaking the dynamic pressure varies from streamline to streain- 
lite, For technical purposes only the average velocity can be utilized. The 
dynamic pressure for the average velocity is slightly lower than the average of 
the dynamic pressures; in other words in defining the dynamic pr -ssure with 
average velocities one takes into accsunt the mixing losses. Thus, we have accord- 


ing to Bernoulli's ecnation 


The efficiency of the process in the pump is best defined by 


eo Pa 


Ea 


n 2 
where Q=2rr., 2duy denotes the amount of fluid entering at T. Figures 29 thru 


i 


31 show, for w= p21, the moment M, the total pressure Pa and the efficiency 


versus ~ for the chosen values of &. and d. 


r 0 
i 
According to the linearized approach the profiles at laxge values of r do not 


depend upon the flow field at smaller radii, thus dynamic pressure and torque are 


determined only by the exit cross section. However the stztic pressure arises by 
15 


an integration over the pressure gradient in the radial direction. The losses 
encountered appear in this analysis as Icssee in static pressure. 
The exit cross section does not play a special role in the computation; in 
other words, each value of r can be considered as the other radius of the rotor. 
In the vicinity of the excrance, one has always a slight pressure drop, for 
large volume flow the pressure drop is more pronounced, and this fact is of 
interest from the point of view of computations. High mass flows are also detri- 


mental from the point of view of efficiency. 


APPENDIX 


For the numerical solution of the system of equations (4.2) we set up a net- 
work of lines x=constant and y =constant with a spacing of h and @ respectively 
along the x and y axes. The length ( is chosen ay une twentieth of the distance 
from the wall to the median plane of the disks, The derivatives in Eqa (4.2) are 


replaced by difference quotients, e.g. 


ba Fite ik ow — “itt, tl ~ “itd, kel 
bk] th ed 


aR h oF 2 

(Al) 
ata | Tien, ett "7 Mien ee * Sie eed 
ay2 


Here i and i+1 correspond to successive lines x = constax’ and k~-1, k, and 
k+1 to success: e lines y = constant with k=0 at the wall. starting with given 
values of i, V, W at che line x= const (i=0) corresponding to the inlet, we 
determine the unknown values of u, Vv, W on the (i+1)st line from the known 
values on the i-th line. 

It is convenient to introduce as new unknowns the differences ai, av, 


aw defined by 


OP Ba oS 


and similarly for Av, and Ow,, and neglect all terms in the equations of higher 


k 


order inthe A's, (Naturally these A's will also depend on i, but only one value 


of i is needed at a time), 


Equations (4,2) then lead to the foilowing linear difference equations in the 


unknown differences Ai, AV, and AW and the unknown pressure function f(x). 


{ - soak = toe L 
(he FF aH ea gl TH te en ee tT Hy! 


20° i Kenav- f(x), 


{ . - = h 
(l-3- wl t2¥,, $m V a (1+ ) (A2) 
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-V arta tee’. Ly yu 


h i, k+l i, kel 


e - 
+ Vs eat Ve aw, + ial & +2)ad, 
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+(2+— u lit ay, + he ya 


“ik k-1’ 
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In addition to Eqs (A2) the complete specificatinn of the problem requires 


Eqs (A3), (A4), and (A5) resulting from boundary conditions. 


Since G = V = W = 0 zlong the wall, one has 


Au, = AV) = Aw, = 0. (A3) 


Since & and V are symmetric and w antisymmetric with respect to the 


median plane of the disks, 


21 19 
4V,,-49,,=0, (A4) 
459 = 0, 


In the continuity equation there appear only first derivatives with respect to 
y, while the other . 0% ‘ons contain second derivatives, C. ll notice that in 
Eq (Al) the firs: derivative of W with respect to y has uot been formed from 


values Ww at adjacent points, This was done in order to obtain a central difference 
2- 2= 


9% 
formula for — » just as for _ and 2x » Asa consequence or.. equation 
by ay dy 


is lost. Therefore, the following approximation is made in the vicinity of the wall. 


an - 
In the last of equations (4. 2) the term ra is zero at the wall; u may be approxi- 
x 


mated by its linear term, Thus 


From the last equation one finds 


and so 


An alternate form is 


Here the second term ought to be zero, but by proceeding in this fashion we are 
more sure tha. : .*: equation of (A5) is not violated bec . of error accumu- 
lation. 

By numerically integrating the radial velocity u with respect to y for any 
fixed value of x, we obtain the total fluid volume passing through the cylindrical 
surface characterized by x, and comparison with input volume provides some 


check of the numerical results. 


Fig 1. Schematic diagram of the flow between two rotating disks 


sutstousesestessesestafostasestastaraataseatastast natestast 
ioneee Se2 c00S0s5R8 aEan : csces! ae — 
my beesoae neue pewecesed teeteceeed | EEEEESEE 

iaaitieen ities 


as 
+ 
rh 


SS 

an. “eRe 

Baaee 
aae 
an 


——p- ad 
versus ad 


Fig. 


2 k, and k, 


2c 


SRRRSRR a ee Lp ae Races 


seepagugeeseeeegeuesseaGeenceusesussere senSUUDeeeusececuas seseuees 
iiteny : Ssiecaereattoaiiiait rH 


an — 


tt : 
pti rttyioththh 
So a ° a roe) .o + 

« e = es s se s 


‘a SRRERanaanavas reer 


STEEEERN 


+L 


{- 


+4 
an 
alo 
ri 


ea 


cas od abe ae Z : ‘ ‘ 


AVAILABLE COPY 


4 


} on 4 
vey: ah 


pair od 


{ 


i pee oe ere ee 


ype 


pasha 


tot 


oe oe ae es ory 


So 
ii Senn AReseaek ames 
PERE By SGUHGSRURGRSUNeGesageeees © IEE ERES 
QuannyEy yee Coe BERSRSRRRenESs b cnes ines 
guauae geeeess Cee So Sanee ¢ (Gaeieome 
PCE oe Ty s CO 
oe OG A a a rer Tr BREE EE HE 
1 SeEyesuens COL Hs RUSEaEE CCC 
gas saannn aannee Le Bean Buel s isme Saeee 
~ Eee tH 4 Ht 4 MBaS!i Fins Ree 
a ee ATE one CA 
nite Sennan: t+ 
1 rT umn 88! a 
sauneen aunGaEes . mepeaal: Too 
. guuuanss AA 
Corr ane CO O) 
CCC corti anti a) Se ane 8 
CCl coat S001 2h. ee 
an aa gaan BER a an 3 
% au Sune am an 
ry naan ree a 
|) Treo ir Baa o 
cco on 3 
cI Seennaenun ae 
Seeuenans an 
Coo Corr 
SUSsnanane Th ® 
seecaeneeee \quaeeus FA 
os ry uae % 
i (82 R00 
3 TORTI g 
‘ Suanunun ‘3 
rtty a 
corey ba 
mme > 
anae 
a oneas es 
+ ann e o 
° aan Beauenae A 
° | Baa 8 
7 a SO re oo 8 
SERGRQUGRY QuEBEy Gua d\\cnonsees ny 
a UORARe Ae ws ee > 
Trt HH CIT} aeneaene RPE EEE HAH = 
Bsasanarescusasbet'd bestasecataesestay atate ssstece itecares aE 
’ fi 
6 SS eg 
3 a Bo Mista! 0) (aly tf s semua, a oy 
’ i cl it 1 . 
tty : aoe 8 rr 7 2G8N8 REE RBNeees | 
qeeee tthe a fe CHEE, aouey duaneae rot os > 
ity] + oes eee : J ] $ 
BAG ANG Re einem i es 3 
Raion mbog eae fare a 
. ies es ) 
= 0 a ed BP a Te bel 
° a a 0 
, Lt +4 44 S 
aeannecuans 45 ; 
rrr er 7 
~ — “tt 
oan ene a eee At Fa 
~ aan \ i ¥ 
oct ee HHH ty mt 4 a} 
Bie ccsecens:! copeagee : 
eeeteaeey jocreetel 
SSeGRbe RO) NANSSRRARY anuy Tt rH ‘4 
tpt fad -f- dott 4 -t' 44d | 
prt tt tte 4 pe RA ee G0 nnn SERSaEE : 
= Date th y et 
2 SEE eg 
q ae: cépecacer>ceegecereysocensts MEE 
' fj 
ze 4 ca e 
; ane, to a 
b A coop rm 
a 
e 
[] 
‘ 
lo] 


negness 
neeeReS: 
seen 


See 


gee tessa tes 
si Periseed 


a 
ann 
ana8 
‘ae 


sea 


® tte 
i . mee 
ama WV 


WI 
pe Tt 
iN Ty 


secerenre 
ATH 


ae 


ae x 
Hic EPPPEEEE 


saisesterisestaaiiit 


dageeceeseesccst 
BOSeaenoounedeans an RON 
auauaeeens 
Sages CR GRSSRaauas Sdugevqnae.uane a 
BRR HS NaSentaaeen tone 
roy ECCLCCE ALLEN GELLER ENG ose 
aw an San \eebeven CRT REET 
2 anne : TT 
ry a8 
annene 
ro a 


San" 4ene.6 


SueGe 40nesaune 
RE SSGG/. Shr ake caee 
‘4 


Peace.s ae” 


tt Aree 


ee 


ea 
an 
a 
‘ste 


-0,2 ~OQ,1 


-0. 3 


-0.5 


-0.8 


-0.9 


-1.0 


zid 


a— z/d 


-0.8 -0,7 -0.6 -0. 5 -0.4 -0.3 -0,2 -0. 1 


-0.9 


+1,0 


See eee 

suggare 
aeee 

as 

wae 


ae 


VOR ER SEES: 
Baa ee 
Se 
bt tat iz 
a 
i i 


aw TGRSe 
NSESRe AREER 


Ize 
gan 


ee iT 
-CBTsane 


v vu 
~ ~ 
N S 
TTT TITEL Le 
Saree 
-_ ~ 
J os 
t r 
“ 
6 
ce) Nn " 
° eo Io 
: ’ 
{~] 
e 
mm 
" 
¢ 
(a (0) i) 
@ 
© rn) 
; 5s 
[> 
- 
3 b: TTT Ith g 
é 3 1 ‘aekees r 
‘ ‘ wyaesn/anas 
ae q 
& 
w w q 
3 é 
] ' 4 
oO we) 3 
° e 
6 o 
’ ’ a 
é w i 
es e 
q 3 3 
ot 
“ 
e ry 3 
e e 
3 
§ 2 
ad r] 
teCWERNaC 
aa a o 
a e 
ro) ° 
8 iy 


-1,0 
-1,0 


rama 


—— a a 


sasacaaas 


. R88 ang—> *FERSRR 
BBe eee sans — SSP SIT TSE aR S 
+ be soe tae om at ab tot oo ot ee on oe ee = SUSE SHOTS CHAED RSE AR Ee aes tase 
i SRE SSSSOMNPTR SAS PHS SSCRKSDS SAPS ESASS SHEEREEAEES CURB AB eenenswaes 2 * TT 


0 T‘O~ 2°O- €°o- ¥“Oo- s°o- 9°0- 2o- 8 °O~ 
P/t 
pike 

0 bt * bad z°o~ € “o> v°O- s°o- 9°O- 2°o°- 8 °0O- 


V 
mnie 
La 


pat--E— ot eT te tt 


\am. (SBR eeaait 


Vu LIS! 


nie 


ett tt 


6°0- 


a 
SRBREBARRE WAR *S 


"o- 


26 


{ 
{ 
' 


-0.1 
-0.1 


ae zfd 
aw————— zf/d 


0.2 
-0.2 


-0.3 


-0.3 


-0.4 
-0.4 


emilsseeen 
eel anne se 
PE Sa 


-0.5 
-0.5 


-0.6 
-0.6 


ae 
eee Ca Weenie 
A OER tateee ae 


-0.7 
-0.7 


Lil HAY 
Ht PAA 


-0.8 
-0.8 


ann 

sane 
-0.9 
-0.9 


annee 

gene 

Po 
-1,0 


Bu 


weeaees 
a 48 
oa 
HH 
eum 
@ awa ue 
adeet ae an “oe 
ganar paar e eapen 
wine sea aseeean 
aa es a pane 
evn 
pice 
ena r 

gaun a 
HH om 
HH oun 

iT yiltlltitt itt 

oa tate Ed EEE 

0 SERRA see eeenes ARaee 


sabaeenes 
2) see eeees dees 
? taht te tt ta 
ret ats ttt 4 
BSGne RRR REEY BS mms 
tA PCE at 
Fh aL dendend gl te 4 


ry 
. 


we 
szecsqsS 
o- 

Ht 


-—— 


: 


i 2sa2f* sae 
anne 


[eee dnd pn 
At 
waue asaeeeeee 

a aranee 


7 


oe wn oe Oe oS > 
ae ae 


@emema 2 
« 


ay 7 


Yo 


ty V for 


: 


da relative tangential veloci 


~-ity 0 an 


2 


Fig.6 Radial ve 


onguasecsagae tai ss 
beaut 


2.0 


eunnecernsea | bended 
quencenaue sit ams 


nat 
20.1 


Gas 
cape te /d 


tee 2d 


<e t 
qunag 


puree 
: 
{ 


oa 


peaad 


sae22e! 


pital : sa uee 
PTL eet ttt 
4 SERS 
‘ yt 
+ - i 
tot 
“0, 2 


eke 
tse 
Lit 
es: 
tt 
SE 
ett | 
ris | | 
Pott {4 
=I eee 
: He 
Hid 
20, 2 


ey 
ra 
aue 
co 
i 


wy eee Sess 


ried t 


neers 
eeu 


eae 
ane 


. 


H 
ne 

as | 

0,4 


i 
au 


Po La PE at 
TLE 
Thi i 
Tile 
sacceatl 
cid 
face? 
seneke 
_ -2808 

tt 
@ 
agen 
-0,5 


rt 
+t 
27 


1,50 
é 
oan 


mena wen 
Ta 
te 


+0, 6 


Pet 
iy 


f 
1,10 


3-- 


2 aape 


Radial velocity u and relative tangential velocity V for ug = 0.25 


= 4 " 

cary 74 on 
SHE 22 WEAN : 
eecinen a — S rom b Ree o 
pee bbe dy “s 
pasesece ote s 


oi 

zee 
eae 
1 


esas 


if 
fe 
HAR AnD AEE: 
f 
? 
0,7 


Fig. 7. 


at: ¢ foe ee 
‘ 7 ay 4 ema ae 
i ee ae 


ESS 
SeseEe: maith: ti 
Sea sem 


0.6 +05 0,4 20,3 “2 0,1 


-0,7 


-0, 8 


aid 


— 


«a——— 3/d 


0,8 =0, / 0,6 0,5 20,4 =u, 3 0,2 -0,1 


20,9 


21,0 


BBS 68ER2 ChSeRs EEK SUEESSE aE 
$ ¢658a 


tig 


iit 


Ni 1 
abe 
. r 3605 


geese 
{ 


Z 
4 
a 


oft rae 


Ty 
HI ae 


+0,2 


0. 


a: 


= 01 


u and relative tangential velocity V for 


~ 


Fig. 8 Radial velocit; 


38 


HAE ERETN PE 
SCT re aut seas fy 2°07 
Pa ia ry 


PH eH —H 


0 t*o- 2°07 € ‘o- +°0- s‘0- 9°0- L°o- 
p/z 
[2 _ 
0 tte z°O- €°O- ¥°O- Ss 9°0O- L°0O- 
BEGGS Reee ee HH FH Pett sen sagas 
SEaun on a rriitisitiriritiT yy Y | gaan sane ‘2 aR 
SEBEB EReue BESSEE ee eeeee sa Zanes ake ‘2 28 
SER Been EES ne SSR 9 8eeR Baseagunesese Py sauna a 
SEREEREREERSASREER SER SERREES BERR see an88 GEsaka 
SER RGRaSEE rT ir y PTT ityiri yr yy ane ae Sa RaBa 
rrtititiiryrr errr yar a augenege a BB saysen 
BEae rtijs Citrtii ri iiiih Ht HH Seeun anus Et ela 
TTT LLLLE EO eee anee a an 
TCT TIL eo oe BaBs aa i 
At Prete +i ti iy trp sane aae8 segue 
REUSSGRREEE Re BSGeSE BRERESESRe ee San Buae 
SRR ERSEERERERREERRSE BERR eee suRe ry ty 
QRS SR ESES SAN ERE EEE eee ee eee aaaeneen cece 
rrTiie”” 6vrTLELLLLL Le “ uss 
sa58 Pitty ity PH Tt Ban aanee ase ‘2 YT 
aaa suneee seeeeseueauue HH rterTy TT ty 
ft saaeeeeee HH rriyiy iy SERB SESERER 
Sanne =f puneesenes Ree asascenanaaaaues 
Pri ti ti iy pee Saugsusesaeuseas 
SEegeeeee Re ee oe SBERERSEBeeSSS 
SEGRSUSEER SRB ESRSES Bese selesseseeanene 
SEReanR Ree COCO ar. + 
SRRREEEEER rivtetiriryy rir yr i VT | 
SUBRESEAES ERRRREE RE RSSRRAAE RASA SSSe— 2 cif mia 
SSERBEEESE 4 {| SSSEERERESESESeK- aceess 
HH HHH BSaREESeReeS=-.4. 88 ae pi 
tosnee fasesecs=-sa,cueeune8 Beuageea maa 
ttt biel —-t—t- $4 -sonSGG00 S0SSSeReee see at rade 4 L, 
seccgeeaeseeecs spit at ate 
ie ‘a 
28 SSSR Peske sae s vi tT 
BopanEeens SERuUSSEEEED: 
-- -—S 885005005555 aes tA 
a2e88 B saner.sen 
secaceees SEeeRe oe 2en 
Ss 2aeReeeRe A ean. 
rt ‘2 HAG 
| wT 
a PEE 
au Seuguaueeee eau 
BB Sas CRASS. La Bee: 
icabescecee an Pi sen 
Baesuusasear aurann rii 
suecssuesessas caaees eeseet 
Segue «aagnene Ssezeaanan*. Hatt eg eases eee ssasens or 
Ae aoe 2 ttt SEenREESs 
FlrASES Saas Bnr. snenriaees ea =GEEEEaR SEnaua 
rey yi dy SuaeP2s S59 2260-4488 Suaegeannn eens: 
ry yt Ba aeeer.eee 225°" .28 08 SERER EERE aSS asgess 
BCL PEER EERE EEE EEE EERE EEE rh 
HH H S22 ascaaae Reet SERESARRRARERREE ee 
222. AES: SESEBeaae a Bagge 
py Sree sae a Pt tH Ht Ht ——— *T 
2 Seung HHH pip ry | 
PT ‘= SESE ERR a8 
a 2enue an [i a 
25 BERS an 2 pty 
sana GRR ee BESaReReneR a e aaen 
SRSERRER FERS $+ PTyTpTr ppp pera pry gy 
a SSSaRk SaSeERAs BEREEESESESE SEnee a 28 aa Be. 
aauase 4+ SERB ESSRESERRRSREE eS aaa aSsEenR gennee a 
TTT L Lire rec aa a aa Suunepaeesee eT 


PTT cEnOSaRSER EE Be 
HH FSSSRSGRRRGt GG Hi 
a ge SESEUSESER EEE is 
ct =z 2 SERUREERESEEE a 
sane uN H H 
asan mi 
Ban rn 
aaae a a 
rit. os “ an COCCI ee B 
non ee 3 COCO T eae s. 
5 g $ Li i 
. t He I 
. : ¥ 
: i] ] t 
i 
ry 
ie t 
2 N ad | 
J 
y ¢ 
Tee r 
ttre r] 
(aie) * ch 
° 
9 gS on 


v ‘ 
J e 
°o i~) 
4 t) 
w w 
e e 
o i) 
8 t 
Byes 7 
wv.) Oo ’ tsa 
J ° K 
3 3 “t: 
1 ' * 
SESE eesaeee Hii 
suaceReaue rr agen .en 
sescebeeaesrast HH i 
i f] 
aa dunes 
Lad ~ f J Pie tt eo eae 
o Pa fa 14 a Tl | Se sennatn seen 
aeRen ee ‘ : 
CCC 
NTEEEERE LLL 
; i 
Sean Se co 
¥ i 9] [- °) 
o o 
¢ 
auar. 
VTA FCT tA 
on a een 'aerrar 
2 28 £9.49 sau ee 
S ° Cyr 
' ' writen 
Vy Ho 
HRS 
a 
° ° rr a 
we Tl dlededed Cy 
' 
—) N wv oO eo °o 
2 a 
° So ° So 
‘ ’ —>|> t t 7 


Fig.9 Radial velocity U and relative tangential velocity V for 
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